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Motivated by lattice QCD studies we investigate the RW transition 
endpoint at imaginary chemical potential in a two-flavor PNJL model. We 
focus on the quark-mass dependence of the endpoint using different forms 
of the Polyakov-loop potential. 
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1. Introduction 

At imaginary chemical potential QCD has an interesting symmetry, 
known as the Roberge- Weiss (RW) symmetry [T]: As a remnant of the 
Z3 symmetry of the pure SU(3) gauge theory, certain shifts in the imag- 
inary chemical potential [i = i9T can be undone by a Z3 transformation 
leading to a periodicity of 9 — > 9 + 2-7rA;/3 with integer k in thermodynamic 
quantities like the pressure. At large temperatures the system undergoes 
a first-order transition jumping between different Z3 sectors when cross- 
ing 9 = {2k + l)vr/3 for fixed temperature. Due to the periodicity #-even 
quantities show a cusp, whereas 0-odd quantities have a jump. For low 
temperatures this transition is a crossover. In between there must be an 
endpoint of the RW transition which can be of first or second order. If the 
transition along 9 = tt/3 ends in a first-order transition, there must be first- 
order lines departing from it implying that the endpoint is a triple point. 
As first-order phase transitions and second-order endpoints might influence 
the phase structure at real chemical potential, this warrants further studies. 

Recent lattice QCD simulations at imaginary chemical potential for two 
and three quark flavors have shown that the order of the RW endpoint 
depends on the quark masses [2 [3]. For low and high masses, the transition 
is of first order. A first-order transition at large quark masses is to be 
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expected from the limit of pure SU(3) gauge theory. In the intermediate 
mass range the transition changes to second order with tricritical points in 
between. 

Since lattice studies are hampered by the sign problem and are com- 
putationally very demanding, it is worth studying these aspects in effec- 
tive models. In the Polyakov-loop extended Nambu-Jona-Lasinio (PNJL) 
model, which can be applied for real as well as at imaginary chemical po- 
tentials, we thus investigate the phase structure in the /x 2 — T-plane. At 
imaginary quark chemical potential the PNJL model also features the RW 
symmetry and we find the RW periodicity as well as the RW phase transi- 
tion. 

The PNJL model at imaginary chemical potential has already been in- 
vestigated by Sakai et. al. [HE]. In a two-flavor PNJL model we extend 
their work and study the order of the RW phase transition endpoint for 
different Polyakov-loop potentials and analyze its dependence on the rela- 
tive strength of the potentials [6j. This is done in two ways: Since quarks 
with larger mass have a smaller contribution to the pressure, increasing 
the quark mass makes the gluonic part more important. Alternatively, we 
directly change the prefactor of the gluonic contribution. 

2. Model 

We employ the PNJL model for two light quark flavors at real and imag- 
inary chemical potential in mean-field approximation following the standard 
procedures. The Lagragian is given by 

with quark fields if), covariant derivative D^, bare quark mass mo, coupling 
constant gs of the four-quark interaction and the Polyakov-loop potential 
U, modelling the gluonic contributions which depends on the Polyakov-loop 
variables <£> and <I>. Parameters are taken from pQ. 
The extended Z3 transformation [5] is given by 

9 + 2irk/3 

$ -> $exp [-i27rfc/3] with k € Z. 

A convenient definition is the modified Polyakov loop, ^ = <3?exp [iO], which 
is then invariant under the extended Z3 transformation. It can easily be 
shown that the PNJL model is invariant under the extended Z3 transforma- 
tion and thus possesses the RW periodicity at imaginary chemical potential. 
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Fig. 1. Modified Polyakov-loop variables and the normalized chiral condensate ct/cto 
at 9 = (left) and 6 = 7r/3 (right) as functions of temperature. The phase of $ 
vanishes at 9 = and only the positive branch is shown at 9 — ir/3. 



3. Results 

We start with a logarithmic form of the Polyakov-loop potential [7], 

= _o(£)qq + 6 ^ log ^ _ + 4 ^3 + |,3) _ 3 ($§)2j _ 

We show the behavior of the order parameters at fixed 9 = and 9 = it/ 3 
in Fig. [1] Along 9 = ir/3 which is in the middle of the period we find a jump 
in the absolute value and the phase of the Polyakov loop. The dependence 
on 9 at fixed temperatures close to the RW transition is displayed in Fig. [2j 
At temperatures higher than the transition temperature Trw the phase has 
a jump and the absolute value a cusp when crossing the RW phase transi- 
tion, signalling the jump from one Z3 sector to another. At temperatures 
slightly lower than Trw we however find two jumps in the phase and also 
in the absolute value. In addition the chiral condensate picks up the same 
discontinuities as the absolute value of the Polyakov loop. For even lower 
temperatures all transitions are continuous. We summarize these findings 
in the PNJL phase diagram shown in Fig. [3l Crossover lines are determined 
by the inflection point of the Polyakov-loop absolute value as a function of 
temperature. Chiral crossover lines are ommitted as they are not relevant 
for our current analysis. 

Using the logarithmic parametrization we find the RW endpoint to be 
a triple point independent of the quark masses, contrary to lattice results. 
First-order lines departing from the triple point are nicely visible. Increasing 
the bare quark masses tuq leads to larger effective quark masses. This results 
in growing "RW legs" , see Fig. [U For mo larger than about 180 MeV the 
first-order lines even reach across the /i = axis. This scenario is shown in 
the right panel of Fig. U]in comparison to the standard-parameter results. 
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Fig. 2. Dependence of the modified Polyakov-loop variables and the normalized 
chiral condensate on 9 for different temperatures around Trw = 190.3 MeV (red 
solid: T = 185 MeV, green dashed: 188 MeV, blue dotted: 191 MeV). 




Fig. 3. Phase diagram in the 9— T (left) plane and the fi 2 — T (right) plane. Red solid 
lines denote first-order RW/deconfinement transitions, green dashed lines show the 
deconfinement crossover, and the blue solid line at real chemical potential denotes 
the chiral first-order transition. The diamonds represent second-order endpoints. 



If instead the coupling constant g$ is increased, which likewise leads to 
larger constituent quark masses, the same effect is found [8]. 

Next we analyze the behavior of other Polyakov-loop potentials. Though 
all parametrizations are designed to reproduce pure-gauge lattice thermo- 
dynamics their effect on the RW endpoint is quite different. 

The polynomial parametrization [9j leads to a second-order transition for 
all examined quark masses. The reason is, that the polynomial parametriza- 
tion shows a much weaker first-order transition in the heavy-quark limit. 

Similarly, the Fukushima-type Polyakov-loop potential [10], given by 



= -bT (54e- a/T $$ + log [l - 63>$ + 4($ 3 + $ 3 ) - 3($$) 2 ]) , 
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Fig. 4. Left panel: "RW legs" in the 9 — T phase diagram for different values of 
the bare quark mass mo- Right panel: Phase diagram in the — T plane for 
two different values of the bare quark mass. Red solid (blue dashed) lines show 
first-order RW/deconfinement (deconfinement crossover) transitions for a high bare 
quark mass mo = 200 MeV. Thin lines show the RW, deconfinement and chiral 
transitions for the standard value of uiq = 5.5 MeV. 
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Fig. 5. Temperature and order of the RW transition as function of parameter b. 



produces a second-order transition for small quark masses and changes to 
first order only for very high quark masses where the PNJL model is not 
applicable any more. An alternative way to drive the system towards the 
pure gauge limit is to increase the global factor b of the Fukushima-type 
Polyakov-loop potential. As presented in Fig. [51 the RW endpoint changes 
from second to first order at about b = 0.09A 3 whereas the default value 
for N f = 2 is b = 0.015A 3 . For b > 0.5A 3 the "RW legs" reach across the 
temperature axis. With increasing b the transition temperature approaches 
the heavy-quark limit of T c = 270 MeV. 
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We conclude, that the PNJL model together with currently available 
parametrizations for the Polyakov-loop potential is not able to reproduce 
the mass dependence found in lattice QCD studies. Sakai et. al. have shown 
that the 'entanglement' PNJL (EPNJL) model, which uses a Polyakov-loop 
dependent coupling gg, reproduces the desired behavior [TT] . 



We have shown that the choice of the Polyakov-loop potential parame- 
trization has an important influence on the order of the RW phase transition 
endpoint. Modifying the strength of the quark degrees of freedom relative 
to the Polyakov-loop potential which models the gluon degrees of freedom, 
we find interesting changes in the phase structure at imaginary and real 
chemical potential. Results from lattice QCD should be used to constrain 
the Polyakov-loop potential parametrizations used in model studies. 
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